基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

通过表面机械加工调控Zr52.5Cu17.9Ni14.6Al10Ti5非晶合金的结构和韧性

作为一种新型结构材料, 非晶态合金的韧性需要进一步提高. 提高非晶态合金韧性的方法有引入枝晶相、调整其成分改变其泊松比影响其剪切带衍生、裂纹扩展等.本文通过表面机械加工的方法来调控非晶态合金的微观结构及韧性. 我们采用真空电弧熔炼、亚稳态薄板离心浇铸系统制备了Zr_52.5Cu_17.9Ni_14.6Al₁₀Ti₅ (原子百分比) (Vit105)非晶合金板,并用表面机械研磨处理方法(surface mechanical attrition treatment, SMAT), 在Vit105板上形成纳米尺度局域类晶体序结构. 基于差示扫描量热分析、纳米压痕实验, 我们发现SMAT处理后的Vit105合金板表面附近弛豫焓更低, 微观结构更加均匀、稳定. 通过显微维氏硬度计测试, 发现SMAT处理后样品的表面附近硬度增大,硬度值分布也更均匀. 通过三点弯断裂实验, 可得到SMAT处理后合金板缺口韧度值从70.7 ± 4.7 MPa·m1/2提高到112.8 ± 3.7 MPa·m1/2. SMAT处理后合金板断裂后, 缺口前端剪切带密度比未处理的更大. Vit105合金板韧性的提高源于SMAT处理对剪切带萌生的促进作用. 该研究表明,表面机械加工可以在非晶态合金中形成局域类晶体有序结构, 影响其结构均匀性, 增大其硬度, 促进剪切带萌生, 提高其韧性. 表面机械加工作为一种新型的改变材料性能的手段, 具有广阔的应用前景.      

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A framework to perform high-order deconvolution for finite-volume method on simplicial meshes

In this article, a new framework to design high-order approximations in the context of node-centered finite volumes on simplicial meshes is proposed. The major novelty of this method is that it relies on very simple and compact differential operators, which is a critical point to achieve good performances in the High-performance computing context. This method is based on deconvolution between nodal and volume-average values, which can be conducted to any order. The interest of the new method is illustrated through three different applications: mesh-to-mesh interpolation, levelset curvature computation, and numerical scheme for convection. Higher order can also be achieved within the present framework by introducing high-rank tensors. Although these tensors feature much symmetries, their manipulation can quickly become an overwhelming task. For this reason and without loss of generality, the present articles are limited to third-order expansion. This method, although tightly connected to the k-exact schemes theory, does not rely on successive corrections: the high-order property is obtained in a single operation, which makes them more attractive in terms of performances.     

The Szegö kernel on a class of non-compact CR manifolds of high codimension

We generalize Nagel’s formula for the Szegö kernel and use it to compute the Szegö kernel on a class of non-compact CR manifolds whose tangent space decomposes into one complex direction and several totally real directions. We also discuss the control metric on these manifolds and relate it to the size of the Szegö kernel.     

Construction Theory of Function on Local Fields

We establish the construction theory of function based upon a local field K p as underlying space. By virture of the concept of pseudo-differential operator, we introduce "fractal calculus"(or, p-type calculus, or, Gibbs-Butzer calculus). Then, show the Jackson direct approximation theorems, Bermstein inverse approximation theorems and the equivalent approximation theorems for compact group D( Kp) and locally compact group K+p(= Kp), so that the foundation of construction theory of function on local fields is established. Moreover, the Jackson type, Bernstein type,and equivalent approximation theorems on the H ?lder-type space Cσ(Kp), σ 0, are proved; then the equivalent approximation theorem on Sobolev-type space Wrσ(Kp),σ≥ 0, 1 ≤ r +∞, is shown.     

一种非正交对称条带束流位置监测器设计

介绍了一种条带束流位置监测器(BPM)的设计与仿真方法。在国家同步辐射实验室"太赫兹近场高通量材料物性测试系统"工程项目中,针对波荡器出口处真空室非正交对称性的问题,设计了矩形真空室和跑道形真空室下的两种非正交对称性条带BPM,并与传统的圆形真空室下条带BPM进行对比。基于边界元法,利用MATLAB软件分别对三种真空室下的条带BPM进行建模和仿真。仿真结果表明:相对于传统的圆形真空室下条带BPM,矩形和跑道形真空室下条带BPM灵敏度提高了30%,阻抗匹配误差相对降低了20%,束流位置拟合误差降低了80%。考虑加工精度,矩形真空室下的条带BPM更适用于该工程。